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The ground state phase transition in the distorted S — 1 
kagome Heisenberg antiferromagnet (KHAF) is studied by 
means of the perturbational calculation and numerical exact 
diagonalization method. For strong \/3 x \/3 lattice distor- 
tion, the hexgonal singlet solid (HSS) ground state of the 
uniform KHAF is destroyed and new singlet state, large HSS 
(LHSS) state, which is globally different from the HSS state 
is realized. The quantum phase transition between these two 
singlet states is analogous to the Haldane-dimer transition in 
the S = 1 antiferromagnetic Heisenberg chain. The presence 
of this transition supports the validity of the HSS picture of 
the ground state of the uniform 5* = 1 KHAF. 
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The kagome Heisenberg antiferromagnet (KHAF) has 
been extensively studied theoretically and experimentally 
because of the interest in the interplay of the strong quan- 
tum fluctuation and the highly frustrated nature of the 
lattice structure. So far, most of the attempts have been 
focused on the ground state and low lying excitations of 
the uniform KHAF. In both 5=1/2 and 5 = 1 |] 
cases, it is known that the ground state is a spin singlet 
state and the magnetic excitation has a finite energy gap. 
In the S = 1/2 case, there are a number of singlet exci- 
tations below the first triplet excitation possiblly down 
to zero energy in the thermodynamic limit On the 
other hand, the singlet excitations also have finite energy 
gaps in the 5 = 1 case. The present author proposed 
the hexagonal singlet solid (HSS) picture for the ground 
state of 5 = 1 KHAF which is analogous to the va- 
lence bond solid (VBS) picture of the ground state of the 
5 = 1 antiferromagnetic Heisenberg chains (AFHC) 
In both HSS and VBS pictures, the 5=1 spins are de- 
composed into symmetrized pairs of two 5=1/2 spins. 
In the HSS state, these 5 = 1/2 spins form 6-spin singlet 
state around each hexagon, while in the VBS state, they 
form 2-spin singlet state on each bond. 

As a real material, Wada and coworkers @-||] 
have investigated the magnetic behavior of m- 
MPYNN-BF4 which can be regarded as the 5=1 kagome 
antiferromagnet. Therefore, if the HSS picture of the 
ground state of 5 = 1 KHAF is verified, this material 
is the first realistic example of the VBS-like state in 2- 
dimensional 5 = 1 magnetic systems. In this context, 
it is quite important to check the validity of the HSS 



picture for the ground state of the 5 = 1 KHAF from 
various points of view. In the present work, we check its 
validity by investigating the stability of the ground state 
against the lattice distortion which destroys the charac- 
teristic magnetic structure of the HSS state. 

Considering the structure of the HSS state, this state 
is expected to be destroyed by the lattice distortion with 
VS X \/3 structure and the phase transition to a new 
ground state with different magnetic structure will take 
place for strong enough distortion. This implies that the 
HSS picture of the 5=1 KHAF is verified if the presence 
of such phase transition is proven. Corresponding check 
is well established for 5 = 1 AFHC, in which the VBS 
state is destabilized by strong dimerization [l0| - p^ and 
this observation elucidated that the Haldane phase is es- 
sentially different from the dimcr phase in 5 = 1 AFHC. 
Actually, we find a new ground state (large HSS state) 
for strongly distorted 5=1 KHAF and a quantum phase 
transition between these two singlet ground states. The 
presence of this phase transition supports the validity of 
the HSS picture for the ground state of the uniform 5 = 1 
KHAF. 

It should be also noted that the material m- 
MPYNN-BF4actually undergoes a structual transforma- 
tion around 128K to the distorted phase with -\/3 x -\/3 
structure [ p^ . This is another motivation of the present 
study, although the strength of the distortion might not 
be strong enough to induce a ground state phase transi- 
tion. 




FIG. 1. V3 X ^/3 distorted ka gome lattice. 
Let us consider the x -\/3 distorted KHAF given by, 
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where Si is the spin operator and X]<i j>ea represent 
the summation over the bonds around the type-a (a = 
A, B or C) hexgons, which are depicted in Fig. |l|. 

For Ja >> Jb = Jc^ the 6 spins around each A-hexgon 
approximately form a 6-spin singlet state. The number 
of remaining spins is iV/3. These 'alive' spins form 
again larger kagome lattice as depicted by the open circles 
in Fig. ^ Even in the strong distortion limit, the signs 
of the effective interactions between 'alive' spins are not 
obvious in general, because there are various paths which 
mediate the effective interaction. We therefore explicitly 
carry out the perturbation calculation up to the second 
order in Jb and Jq to obtain the effective Hamiltonian 
for the 'alive' spins as. 
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where Si is the spin operator and the summations , 

<nnB> 

are taken over the near- 

<nnC> <nnn> <nnnn> 

est neighbour pairs within B-hexagons, those within C 
hexagons, next nearest neighbour pairs and next next 
nearest neighbour pairs of 'alive' sites, respectively, as 
depicted in Fig. ^. 




FIG. 2. The strongly distorted KHAF. The open circles 
represent the 'alive' spins. The strongest effective bonds for 
case Ja » Jb ~ Jc are shown by thick lines. 
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FIG. 3. The effective exchange coupling constants between 
'alive' spins (open circles). The filled circles represent the 
'dead' spins connected by the A-bonds. 

Using the numerically obtained eigenvalues and eigen- 
states of the 6-spin hexagon cluster, the effective ex- 
change couplings are calculated as. 



JnnB = 1.190573JbJc 
Jnnc = 1.190573JbJc 

Jnnn - 1.112947 JbJc 

•Amnn — 1.035321JbJc 
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0.519357Jg, 
0.669519J^, 
0.519357(Jb + Jc), 
0.519357( J| 4- Jg), (4) 



for 5 = 1 up to the second order in Jb and Jc. Here we 
have set Ja = 1. 

For comparison, we have also calculated the case of 
S = 1/2 as. 



JnnB — 0.388858JbJc 
J„„c = 0.388858JbJc 
Jnnn = 0.350555JbJc 

Jnnnn = 0.312253JbJc 
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a, we have 
0.057519a2, 



If we set Ja = 1 and Jb 

JnnB = JnnC = -0.009608a2, Jnnn 

J„„„„ = 0.019216a2 for S = 1/2 and JnnB = JnnC = 
0.001697a^ Jnnn = 0.074233a^ Jnnnn = -0.003394^2 
for = 1. In this case, the strongest effective interaction 
is the next nearest neighbour interaction. If we neglect 
other interactions, the whole lattice of 'alive' spins are 
decomposed into three equivalent sublattices of kagome 
type which are depicted by thick solid lines, thick dot- 
ted lines and thick broken lines in Fig. |[ Therefore the 
ground state is expected to be the spin singlet state for 
a << 1. For = 1, this state is the HSS state on large 
hexagons which are three times larger than the original 
kagome lattice if the HSS picture is vaild for the uniform 
5' = 1 KHAF. We will call this state as 'large hexagonal 
singlet solid' (LHSS) state. It should be noted in this 
state that the 5=1 spins around the A-hexagons are 
not decomposed into two S" = 1/2 spins but they partic- 
ipate the six spin singlet state as a whole. In this sense, 
this state is an analog of the dimer phase of the S — 1 
dimerized AFHC in which all S" = 1 spins are not decom- 
posed into S" = 1/2 spins but paired into local singlet 
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state as a whole. In contrast to the latter, however, the 
'alive' spins in the LHSS state are again decomposed into 
two 5=1/2 spins and form 6-spin singlet state on en- 
larged sublattices. Thus the structure of the LHSS state 
is globally different from that of the HSS state of the 
undistorted kagome lattice. Therefore, we may expect a 
ground state phase transition at a critical value of a from 
the HSS state to the LHSS state. 





N=18 



FIG. 5. The a-dependence of the singlet triplet gap of 
distorted KHAF with Ja = 1, Jb = >/c = « for (a) 
S = 1, iV = 18 and (b) 5" = 1/2, iV = 18, 27. 




N=27 



FIG. 4. Clusters used for numerical diagonalization. 




To confirm the presence of phase transition between 
the HSS and LHSS states for — Jq — uJa, the nu- 
merical diagonalization calculation is carried out for fi- 
nite size clusters. The clusters used for the calculation 
are shown in Fig. ^. The singlet-triplet gaps are plotted 
in Fig. ^(a) against a for iV = ISjS" = 1. It is evident 
that the gap has a minimum around a ~ 0.5 which is 
indicative of the phase transition in the thermodynamic 
limit. On the other hand, the physical picture of the 
ground state of the 5 = 1/2 KHAF is still unclear. In 
this case, the singlet triplet gap is a monotonous function 
of a for 5 = 1/2 as shown in Fig. |5|(b) for = 18 and 
27 suggesting the absence of phase transition. Therefore, 
the singlet ground state of strongly distorted S = 1/2 
KHAF with Jq = Jq is continuously connected with the 
ground state of the uniform 5=1/2 KHAF. The pres- 
ence of the HSS-LHSS phase transition for the 5=1 case 
and its absence for the 5=1/2 case support the validity 
of the HSS picture for the ground state of undistorted 
5 = 1 KHAF. 

In general, the expressions of the effective coupling 
constants and (||) show that there is a strong com- 
petition between the contributions from different inter- 
action paths. Therefore, we may expect a rich vari- 
ety of ground state phases for the distorted KHAF. For 
example, if we set Ja = 1, Jb = and Jc = 0, 
we have J„nB = -0.251948^2, J„„c = -CUGSlSa^, 
J„„„ = -0.146518^2 and J„„nn = -G.146518a for 5 = 
1/2 and JnnB = -0.669519a2, J„„c = -0.519357a2, 
Jn„„ = -0.519357^2 and J„„„„ = -G.519357a2. There- 
fore all effective bonds are ferromagnetic and the ground 
state is ferrimagnetic with 1/3 of full magnetization. It 
is numerically verified that the ground state is always the 
ferrimagnetic state with the same value of magnetization 
for Ja = 1 , Jb = and Jc = with < a < 1 by exact 
diagonalization for clusters with A^ = 18. The determi- 
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nation of the ground state phase diagram in whole pa- 
rameter plane is an interesting issue although it requires 
a huge computational effort. 

In summary, we have found that the ground state of 
KHAF remains singlet even in the strong distortion limit 
in the case Ja > Jb ~ Jc- For S ~ 1 KHAF, this singlet 
ground state is the LHSS state which is globally differ- 
ent from the HSS ground state of the undistorted S — 1 
KHAF. The HSS-LHSS phase transition is shown to take 
place at an intermediate value of distortion. This transi- 
tion corresponds to the Haldane-dimer phase transition 
in the S" = 1 AFHC The presence of this phase 

transition supports the validity of the HSS picture for the 
undistorted S = 1 KHAF which is analogous to the VBS 
picture of the Haldane phase [^,|| . 

In this paper, we have discussed the ground state phase 
transition in detail only for the case J a > Jb = Jc, 
because this is the simplest case in which the ground state 
remains singlet even in the strong distortion limit. In 
general, there remains the possibility of many other kinds 
of ground states. The detailed analysis of the ground 
state phase diagram of this model in the whole parameter 
plane will be reported elsewhere. 

The numerical calculation is performed using the HI- 
TAC SR8000 at the Supercomputer Center, Institute for 
Solid State Physics, the University of Tokyo and the 
HITAC SR8000 at the Information Processing Center, 
Saitama University. The numerical diagonalization pro- 
gram is based on the TITPACK ver.2 coded by H. Nishi- 
mori and KOBEPACK/1 coded by T. Tonegawa, M. 
Kaburagi and T. Nishino. This work is supported by the 
Grant-in- Aid for Scientific Research from the Ministry of 
Education, Science, Sports and Culture. 
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